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Abstract 

Several recent methods used to analyze asymptotic stability of delay- differential 
equations (DDEs) involve determining the eigenvalues of a matrix, a matrix pencil 
or a matrix polynomial constructed by Kronecker products. Despite some similarities 
between the different types of these so-called matrix pencil methods, the general ideas 
used as well as the proofs differ considerably. Moreover, the available theory hardly 
reveals the relations between the different methods. 

In this work, a different derivation of various matrix pencil methods is presented 
using a unifying framework of a new type of eigenvalue problem: the polynomial 
two-parameter eigenvalue problem., of which the quadratic two-parameter eigenvalue 
problem is a special case. This framework makes it possible to establish relations 
between various seemingly different methods and provides further insight in the 
theory of matrix pencil methods. 

We also recognize a few new matrix pencil variants to determine DDE stability. 
Finally, the recognition of the new types of eigenvalue problem opens a door to 
efficient computation of DDE stability. 
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1 Introduction 



Mathematical models consisting of delay- differential equations (DDEs), in the 
simplest form 

x{t) = Aox{t) + Aixit-T), Ao,AiGC"^", t>0, (1) 

occur naturally in a wide variety of fields related to applied mathematics, 
such as engineering, control theory, biology, traffic modeling, neural networks, 
mechanics and electronic circuits. For most applications, it is desired that 
— >■ as t — > oo for any bounded initial condition. This property is re- 
ferred to as asymptotic stability. If a DDE is not asymptotically stable, the 
linear model typically breaks down or the modeled system has unwanted prop- 
erties, such as oscillations or the energy content of the system is unbounded 
in time causing the modeled physical object to break or at least turn ineffi- 
cient. Clearly, asymptotic stability is important in practice and numerical and 
analytical tools to analyze asymptotic stability of DDEs is a popular topic of 
research. For instance, large parts of several monographs, which are standard 
references in the field of DDEs, deal with stability of DDEs, e.g., the books of 
Bellman [2], Niculescu [25], Michiels and Niculescu [22] and Gu, Kharitonov, 
and Chen [5]. See also the survey papers |27], [TU], and [IB]. This paper is 
concerned with asymptotic stability for ([T]) as well as more general DDEs, 
e.g., DDEs with delays in the derivative (called neutral DDEs) and DDEs 
with multiple delays. These more general DDEs will be considered later in 
this paper. We mention that DDEs that are not of neutral type are also called 
retarded. 

AsjTiiptotic stability is often described using the solutions of the characteristic 
equation associated with ([T]): 

det(-AJ + Ao + Aie~^^) = 0, 

of which the solutions A are called eigenvalues; the set of eigenvalues is called 
the spectrum. The DDE ([1]) is asymptotically stable if the spectrum is con- 
tained in the open left half plane (see, e.g., [221 Prop. 1.6]). 

We will also consider more general classes of DDEs in this paper. For some of 
these DDEs, in particular neutral DDEs, it is not sufficient that all eigenvalues 
have negative real parts to ensure asymptotic stability. A neutral DDE is 
asymptotically stable if and only if the supremum of the real part of the 
spectrum is negative (see, e.g., p2, Prop. 1.20]). 

Because of these properties, explicit conditions such that there is a purely 
imaginary eigenvalue can be very useful in a stability analysis. In this paper 
we will study explicit conditions on the delay r such that there is at least one 
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purely imaginary eigenvalue. In the literature, there are several approaches 
to characterize these values of r, sometimes called critical delays, switching 
delays, crossing delays, or kernel and offspring curves; see [15, Remark 3.1] for 
some comments on terminology. 

One approach to determine critical delays is to consider the eigenvalues of cer- 
tain matrices or matrix pencils constructed by Kronecker products. Methods 
of this type are presented by Chen, Gu, and Nett [4| (see also [9i Thm. 2.13]), 
Louisell [in], Niculescu [21| (see also [221 Proposotion 4.5]), and Fu, Niculescu, 
and Chen [g||2B] . The works [H] (see also [ISl Chapter 3]), [13] and [B] also 
use a formulation of eigenvalue problems containing Kronecker products. Even 
though these popular methods have some characteristics in common, the ideas 
used in the derivations differ. For instance, Louisell [l9j derives a result for 
neutral DDEs by considering a linear ODE which is proven to share imagi- 
nary eigenvalues with the DDE whereas Chen, Gu, and Nett [5] and several 
other authors depart from the characteristic equation and exploit the fact that 
the eigenvalues of Kronecker products are products of the eigenvalues of the 
individual factor matrices. 

The goals of this paper are: 

(a) to introduce a new type of eigenvalue problem, the polynomial two-parameter 
eigenvalue problem, with the quadratic two-parameter eigenvalue problem as 
important special case; 

(b) to show the relevance of this problem to determine critical delays for various 
types of DDEs; 

(c) to provide alternative derivations of existing matrix pencil methods using 
the context of polynomial two-parameter eigenvalue problems; 

(d) to hereby provide a new unifying framework for the determination of critical 
delays; 

(e) and, finally, to recognize a few new variants of known matrix pencil methods. 

For given matrices Ai,Bi,Ci G C"^", i = 1,2, the (linear) two-parameter 
eigenvalue problem is concerned with finding A, G C and x,y G C"\{0} such 
that 

(2) 

= XB2y +fiC2y. 

V 

There is a close connection between linear two-parameter problems and two 
coupled generalized eigenvalue problems involving Kronecker products; see [T] 
and Section [3] for further details. 

In this paper, we will consider polynomial two-parameter eigenvalue problems 
and show that there are associated (one-parameter) quadratic eigenvalue prob- 
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lems which are very relevant for critical delays of DDEs. Note that the use of 
multivariate polynomials, which are closely related to multiparameter eigen- 
value problems, is not new in the field of stability of DDEs. Multivariate 
polynomials are used in, e.g., [T7], [12], and [TT] with apphcations in [5]; see 
also the summaries in the standard references [25l Section 4.1.2] and [9, Sec- 
tion 4.6]. In this work, we discuss a new natural way to interpret matrix pencil 
methods in the context of two-parameter eigenvalue problems. 

The results of this work are ordered by increasing generality. The idea of 
an alternative interpretation of matrix pencil methods using polynomial two- 
parameter eigenvalue problems is first illustrated in Section [21 In Section Owe 
give connections between certain polynomial two-parameter eigenvalue prob- 
lems and associated quadratic and polynomial eigenvalue problem. These links 
are used to derive the polynomial eigenvalue problem occurring in matrix pen- 
cil methods for more general types of DDEs in Section HI After stating some 
new variants of matrix pencil methods in Section \5\ we end with some conclu- 
sions and an outlook in Section O 



2 DDEs with a single delay 

An important aspect of this work is a further understanding of matrix pencil 
methods. The derivation for the most general type of DDE is somewhat tech- 
nical and contain expressions difficult to interpret by inspection. Therefore, 
to ease the presentation, it is worthwhile to first illustrate the general ideas of 
the theory by considering retarded DDEs with a single delay. 

In this section we derive a polynomial two-parameter eigenvalue problem cor- 
responding to purely imaginary eigenvalues of a DDE, and apply a result that 
will be proved in Theorem [3] in Section [3] to identify that the eigenvalue prob- 
lems are the ones that occur in the matrix pencil methods proposed in [3] and 

m- 

First, we introduce the following (usual) notations: cr{A) and cr(A, B) denote 
the spectrum of a matrix A and matrix pencil (^4, B), respectively; / denotes 
the identity matrix, the Kronecker product and © the Kronecker sum (i.e., 
A®B = A^I + I(S)B). If a DDE is stable for r = 0, then denotes the 
delay margin, i.e., the smallest delay r for which the DDE is no longer stable. 

Consider the DDE 

Box{t) = Aox{t) + Aix{t- t), (3) 

where Ao,Ai,Bq G C"^". We will rederive the eigenvalue problems that arise 
in the following two matrix pencil results. We hereby note that matrix pencil 
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methods are generally stated in various degrees of generality, for various types 
of DDEs. Theorem [1] below is for the slightly different setting of neutral DDEs 
and Bq = I. Theorem [2] as stated here is a restriction of the original result |1] 
to single delays. We postpone the discussion of the more general result in [1] 
to Section 4. 

Theorem 1 (Louisell \W, Theorem 3.1]) Let Ao,Ai,Bi G M"""". Then 
all purely imaginary eigenvalues of the neutral DDE 

x{t) + Bix{t - r) = AQx{t) + Aix{t - r) (4) 

are zeros of 

det((A/ - Aq) ® (A/ + Aq) - {XBi - Ai) ® {XBi + A,)) = 0. (5) 

Theorem 2 (Chen, Gu, and Nett; special case of [4j, Theorem 3.1]) 

Suppose ([T]) is stable for r = 0. Define 



U :-- 



I 

Ai ® / 



and V 







-I^Al -Ao®Al 



(6) 



// the delay margin r* is finite and nonzero, then = min^ — where at G 
[0,27r], LOk > 0, and e"^"^' 
Aie-'">'). 



G cr{V, U) satisfies the relation iuj. G o-{Aq + 



Theorem [2] gives a formula for the delay margin in terms of the solutions of the 
generalized eigenvalue problem involving the pencil {V,U), which represents 
a linearization of the quadratic eigenvalue problem (QEP) 



[l2\Ai^I) + fi{Ao 



+ I®Al)v = Q 



(7) 



2 

for G C and nonzero f G C" . An exhaustive characterization of possible lin- 
earizations was recently given in [2D] • Observe that the matrix polynomial ^ 
in Theorem [1] also represents a quadratic eigenvalue problem; a linearization 
was given in pjj as well. (A linearization adapted to the quadratic eigenvalue 
problem in the matrix pencil method in [H] was recently proposed in [?].) 



Both of the matrix pencil methods (Theorem [T] and Theorem [2]) involve 
quadratic eigenvalue problems ([5]) and ([7]). However, for instance from the 
original proofs of these results, there is no obvious relation between these 
two approaches. We will develop a framework that can derive both quadratic 
eigenvalue problems in a unifying manner which gains further insight in the 
relations between the methods. 



Consider the eigenvalue problem associated with 
XBqx = {Aq + Aie~^^) X, 



(8) 
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for nonzero a; G C". We are interested in the case where there is a purely 
imaginary eigenvalue, say A = iu. We denote /x = e"^"^. Under the assumption 
that the eigenvalue is imaginary, i.e., A = iu, we have A = —A and /I = 
This yields 

-XBoy = (Ao + fi-'A,)y, (9) 
where y = x. Hence, multiplying ([9]) by /i and rearranging the terms, we have 

AqX = XBnX —liAiX, 

_ _ (10) 
Aiy = -X^xBoy -fiAoy. 

Now first, for given Ai, Bi, Ci, Di, Ei, and Fi G C"^", i = 1,2, consider the 
following quadratic two-parameter eigenvalue problem 

[Ai + XBi + fiCi + X^Di + X^Ei + /i^Fi) x = 0, 
(^2 + XB2 + llC2 + X^D2 + A/x^2 + Ai^Fs) y = 0, 

where the assignment is to compute one or more tuples (A, ^, x, y) with nonzero 
x and y. As for the linear two-parameter eigenvalue problem we will call (A, /i) 
an eigenvalues and x ® y an eigenvector. We see that (ITU]) is special case of 
this general quadratic two-parameter eigenvalue problem with just one non- 
linear term and one additional vanishing matrix. As an implication of Theo- 
rem [3] in the next section, we will show that the following two (one-parameter) 
quadratic eigenvalue problems are associated with ffTOj) : 



X^{Bq O Bq) + X{Bq ®Ao-Aq® Bq) + {Ai O - Ao ® Ao)\ {x®y) = 

(11) 

and 

fi^{Ai ® + /i (Ao ® 5o + 5o ® Ao) + {Bo ® A))] {x0y) = 0. (12) 



Using these QEPs, we can now rederive Theorems [T] and Theorem [2] as follows. 
Although Theorem [1] applies to the wider class of neutral DDEs, we can restrict 
it to the class of retarded DDEs by setting Bi = 0. Then taking Bq = I in the 
quadratic eigenvalue problem flTT]) exactly renders (jS]) in Theorem [T], under 
the assumption that the matrices are real. 

Similarly, (fT2|) corresponds to the quadratic eigenvalue problem ([7]) in Theo- 
rem [21 note that (fT2l) gives ([7]) if we replace conjugation with the conjugate 
transpose as follows. Instead of IQ as the conjugate of ([8]), we can also take 
the conjugate transpose of ([!])• (In Section O we will exploit similar techniques 
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to derive new matrix pencil methods.) The resulting quadratic two-parameter 
eigenvalue problem is 



Aqx = XBqx —fiAix, ^^^^ 
^ Aly = -XfiB*y -fiA^y, 

where now y is the left eigenvector of ([8]). The second equation in (fT3|) is the 
transpose of (ITU]) . Application of Theorem [3] to (IT^ yields ([7]). 

In the next section we will prove the theorem that we need for flTTl) and f[T21) 
and also more general results. 



3 Quadratic and polynomial two-parameter eigenproblems and as- 
sociated quadratic and polynomial one-parameter eigenproblems 



First, we will review some facts for the (linear) two-parameter eigenvalue prob- 
lem ([2]), see also [Ij. Define the matrix determinants 

Ao = Si ® C2 - Ci ® B2, 
Ai = Ai0C2- C10A2, 
A2 = Bi0 A2- Ai® B2, 

where Aj G C"^^"^, i = 0, 1,2. Associated with ([2]) are two (decoupled) gen- 
eralized eigenvalue problems (GEPs) 

Aiz = XAqz, 

(14) 

A2Z = ^AqZ, 

where z = x®y. (In fact, these GEPs are equivalent with ([2]) if Aq is nonsingu- 
lar; see [1]). As there are two generalized eigenvalue problems which correspond 
to the linear two-parameter eigenvalue problem ([2]), we will see that there are 
two quadratic (one-parameter) eigenvalue problems which correspond to the 
quadratic two-parameter eigenvalue problem. We will see in the derivations 
of the matrix pencil methods that some methods correspond to one form and 
some to the other. 

To be able to handle two wider classes of DDEs, we will prove two results that 
deal with generalizations of problem (fTOl) : we will make use of these theorems 
in the derivation of matrix pencil methods in Section HI 
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The first generalization of (fTOi) . which we will use for neutral DDEs, involves 
an additional cross term A/i: 



Aix = XBix +fiCix +XfiDix, 
A2y = XB2y +fiC2y +XfiD2y. 



(15) 



Theorem 3 // (A, yu) is a solution of ( IT^ with corresponding eigenvector 
(x, y) then: 

1. X is an eigenvalue with corresponding eigenvector x ®y of the QEP 

X^{Di ®B2-Bi® D2) + A(Ai ®D2- Di® A2+ 

- Si ® C2 + Ci ® B2) + {Ai ® C2 - Ci ® A2)] {x®y) = 0. (16) 

2. /i is an eigenvalue with corresponding eigenvector x®y of the QEP 

l?{Di 0C2- Ci® D2) + fi (Ai ®D2-Di® A2+ 

- Ci® B2 + Bi® C2) + {Ai 0B2- Bi® A2)] {x®y) = 0. (17) 



Proof. We show the first implication; the second follows by switching the 
roles of A and /i; Bi and -B2; and Ci and C2. 

Equation (IT^ holds because 

X^{Di ®B2-Bi® D2){x y) 

= X{Di® {A2 - ^^2 - A^D2) - {Ai - fiCi - XfiDi) ® D2)ix y) 

= X{Di® {A2 - fiC2) - (Ai - fiCi) (g) D2){x (g) y) 

= (A (Di 0A2-A10 D2) + A/i (Ci 0D2- Di0 C2)){x (g) y) 

= {X{Di0A2-Ai0D2) + iCi0{A2-XB2)-iAi-XBi)0C2)ix(g)y), 

where we used that 

Xfi (Ci ®D2-Di® C2){x y) 

= (Ci {A2 - XB2 - fiC2) - iA^ - ASi - fiCi) C2){x y) 
= (Ci (A2 - AS2) - (Ai - ASi) C2)ix y). 



□ 



To derive the matrix pencil methods for DDEs with multiple commensurate 
delays (that is, multiple delays that are all integer multiples of of some delay 
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value r, i.e., = ruk where Uk G N) we will need the following more general 
polynomial two-parameter eigenvalue problem 



k=0 k=l 



m 



k=0 k=l 



in the next section. 

We now prove the following result that will prove useful. Associated with 
this polynomial two-parameter eigenvalue problem is the following polynomial 
eigenvalue problem (PEP) for fi. 

Theorem 4 // (A, fi) is an eigenvalue of (fT8l) with eigenvector {x, y) then 



m 



k=l 



+ J2 ® C2,i - Ci,fc ® S2,,)] (x ® 2/) = 0. 

k=l,i=l 



Proof. One may check that 

Aix = XBix+ fiCix + XfiDix, 
A2y = XB2y+ nC2y + X^D2y 



if we let Bi = Bi^o, Di = ELi/'^^i.^ and d = J:T=i f^'^'^d^k for t = 1,2. 
Application of Theorem [3] yields that 
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k=l 



= [{Ai ® ^2,0 - Bi^o ® A) 

m m 

fc=l fc=l 

m m 

- J2 /"'C^l.fc ® ^2,0 + 5i,o ® E /~'C'2 
fc=l fc=l 

m m m 

k=l k=l k=l 

= {Ai B2,o - Bifl ® A2) 

m m 

+ ® E /^2,fc - E /'^l>fc ® ^2 

A;=l fc=l 
m m 

- E Ai^'C^i^fc ® B2,o + 5i,o ® E f^^'^^,k 
k=l k=l 

m m m 

+( E /^i.fe ® E f'^^2,k - E /^i>fe ® E /^2,fe)] (a^ ® y) 
fc=i fc=i 

which completes the proof. 



k=i 



k=l 



□ 



4 Generalizations for neutral systems and multiple delays 



In Section [2] we used the quadratic two-parameter eigenvalue problem (|TOl) 
to derive the quadratic eigenvalue problems in Theorem [1] and Theorem [2j 
However, in Section [2] we limited ourselves to the setting of a single delay 
DDE. The original formulations of Theorem [1] [19j and Theorem [2] [3], were 
stated for more general types of DDEs, which we will study in this section. 
In particular, we discuss neutral systems in Section 14.11 and the DDEs with 
multiple commensurate DDEs in Section 14. 2[ 



4.1 Neutral DDEs 



Consider the neutral DDE 

Box{t) + Bix{t - r) = Aox{t) + Aix{t - r). 



where Aq, Ai, Bq, Bi G C"^". The generality of Theorem [3] allows us to derive 
the matrix pencil methods for this DDE in a similar way as in Section [2l With 
X = iu and fi = e~*™ we note that the corresponding eigenvalue problem and 
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its complex conjugate can be expressed as 



Aqx = XBqx +XnBix —^Aix, 
Aiy = -XBiy -XfiBoy -Mo?/- 

After applying Theorem [3] we derive that 



(19) 



{-Ao 0Ao + Ai0 Ai) + X{-Ao Bq- Bi0 Ai 

+ Bo 0Ao + Ai Bi) + A2(-Si 0Bi + Bo0 Bq)] (x y) = 0, 

and after rearranging the terms we get 

{XBo - Ao) {XBo + Ao) - {XB^ - A^) {XBi + A^)) {x0y) = 0. (20) 

This is a slight generalization of the eigenvalue problem presented by Louisell 
|19j . since in [19] it is assumed that Bo = I and that the matrices are real. 
Louisell, motivated by a connection with a certain differential equation of 
which all purely imaginary eigenvalues coincide with purely imaginary eigen- 
values of the DDE, suggests that (jl]) can be determined by solutions of the 
generalized eigenvalue problem 



A 



1 01 Bi0l 
l0Bi I 01 



w 



Ao0l Ai0l 
-l0Ai -l0Ao 



w. 



However, we note that this is just one possible hnearization of (jlj); any of the 
linearizations in pTpU] might be considered. Moreover, there also exist numer- 
ical methods for quadratic eigenvalue problems that try to avoid linearization; 
see EHI for an overview. 



The second resulting quadratic eigenvalue problem from applying Theorem [3] 
to (HI]) reads 



l?{Bi Ao + Ai0 Bo) + ^ {Ao 0Bo + Bi0 Ai+ 

+ Ai0Bi + Bo0 ~Ao) + {Ao 0Bi + Bo0 ~Ai) \ {x0y) = Q 



At this point we note that we can interchange all left and right operators in the 
Kronecker products to get a special case of the result in ^13j (where multiple 
delays are considered). 

To determine a relation with a result by Fu, Niculescu, and Chen ^ we note 
that similarly to the derivation of (fT9|) using x and y = x, we can also derive a 
quadratic two-parameter eigenvalue problem involving x and its corresponding 
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left eigenvector y: 

Aqx = XBqx +A/i-BiX —^Aix, ^^^^ 
Aly = -XBly -XfiB*y -^iA*y. 

Application of Theorem [3] yields 

{Ao Bl + BQ®Al)+fi {Ao ® B*Q + Bi® A*+ 

+ Ai® Bl + So ® Al) + i?{Bi ®Al + Ai® fi*)] {x®y) = 0. 

This is a special case of the method in [8j which applies to DDEs with multiple 
commensurate delays, which we be the topic of the next subsection) ^ . 



4^.2 Multiple delays 

We now consider the case of DDEs with multiple, say m > 1, delays. In the 
literature there are essentially two ways to handle this situation. Either the 
curves or surfaces corresponding to the critical delays are parameterized using 
m — 1 free variables, as is done in for instance [T^lGfTS] . In other approaches, 
e.g., [HE], it assumed that the delays are commensurate. 

Here, we will focus on the case of multiple commensurate delays as the pa- 
rameterizations do not yield stability information after the solution of one 
eigenvalue problem. 

Consider the DDE with commensurate delays 

m 

Box{t) = J2Akx{t-Tk). 
The associated eigenvalue problem is 
e-"^Mfe - ASo) v = 0. 

\k=0 J 

As in the previous section we substitute X = iu and /x = e~*™ and consider 
the complex conjugate of the eigenvalue problem. After rearrangement of the 



Note that Bi is defined with an opposite sign in 
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terms and sums we have 



(22) 

Aqv = XBqv - J2 tJ'''AkV. 

k=i 

This is of the same form as the polynomial two-parameter eigenvalue problem 
in ([IH]) with Ai = -A^, = Bq, Bi^k = 0, k = 0, . . .m - 1, Ci,fc = Am-k, 
k = 1, . . . ,m, A2 = Aq, ^2,0 = Bq, B2,k = 0, k = 1, . . .m, C2,k = -Ak, 
k = 1, . . . ,m. 

Theorem H] and several manipulations of the sums yield 



-A^(g)Bo + /i"(-5o ®Ao) + J2 IJ-\-Am-k ® Bo) 

k=l 

m 

+ E/^™^'(-^o® A)] {v0u) 

i=l 

m m 

[ - y: /"""'(A ® Bo) - 5:/i™+^(i?o ® A,)] (v ® u), 



k=0 i=0 



which is the eigenvalue problem in derived in [T 

As in Section [2] and the neutral case in the previous subsection, we may con- 
sider the conjugate transpose instead of the transpose. This yields the poly- 
nomial eigenvalue problem in [U Thm. 3.1]. 

Finally, the most general result is for neutral commensurate DDEs. We show 
that the eigenvalue problem in Fu, Niculescu, and Chen [8] also is a polyno- 
mial eigenvalue problem that is connected with a polynomial two-parameter 
eigenvalue problem. Although the analysis is similar as for the previous cases, 
this general case involves more technicalities and more involved expressions. 
Consider the polynomial two-parameter eigenvalue problem corresponding to 
the neutral commensurate DDE 

m m 

Bkx{t -rk) = Y ^kx{t - rk), 

k=0 k=0 



I.e. 



Ao-Y >^lJ'Bk + E /A N = 0- (23) 

^ fc=0 k=l ) 
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The complex conjugate transpose is 

Cm m \ 

/i-A* + Y: Xl^-'-'Bl + Y: r-'Al ]u = 0. (24) 
k=0 k=l / 

We can now combine fl23p and (1241) into a polynomial two-parameter eigenvalue 
problem 



m 



r 'm (25) 

,k A* 



m—k 



U. 



k=0 k=l 



This corresponds to ffTSl) with Ai = Aq, Bi k = B^, k = 0, . . .m, Ci^k = —Ak, 
k = l,...,m, A2 = -A*^, B2,k = B^_i^, k = l,...m, C2,k = A*^_f^, k = 
1, . . . , m. Theorem m yields 



{Ao ® B*^_^ + 5o ® A*J 

m 

+ J2 /«'(Ao ® + Bk0A*^ + Ak® B*^_^ + So ® A*^_i^) 

k=l 

m 

+ Yl f^'''-\Bk0A*^_^^ + A,®B*^_,)]{v0u) = O. 

k=l,i=l 

We note that with some effort it can be verified that the matrix coefficients Qk 
in [SI Thm. 2] are exactly the matrix coefficients that occur in this polynomial 
eigenvalue problem. 

5 New variants of matrix pencil methods 



In this section, we introduce some new matrix pencil methods, which are 
variants of existing approaches. For ease of presentation, we will state the 
results for neutral DDEs with one delay, but all methods can be generalized 
for DDEs with multiple commensurate delays. 

Moreover, we will only mention the relevant quadratic two-parameter eigen- 
value problems (which will be polynomial two-parameter eigenvalue problems 
for DDEs with multiple commensurate delays); as we have seen before, every 
such two-parameter eigenvalue problem has two associated (one-parameter) 
eigenvalue problems, one for A and one for /i giving two possible resulting 
matrix pencil methods. (For DDEs with multiple commensurate delays there 
is just one associated polynomial eigenproblem, for /u.) 
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The quadratic two-parameter eigenvalue problem for the the neutral single- 
delay DDE given in (fT9|) is just one of several possible quadratic two-parameter 
eigenvalue problems. We can get the following expressions by transposing none, 
one, or both equations: 

a) (HHD; 

b) f|T9|) but with the first equation transposed: 

AqX = XBqX +XiJ,Bfx —fiAjx, 
Aiy = -XBiy -X^Boy -Mo2/; 

c) (fT9|) but with the second equation transposed: 

AqX = XBqX +XflBiX —flAiX, ^^^^ 

Aly = -XBly -XfiB^y -fiA^y; 

d) and (fT9l) but with both equations transposed: 

AqX = XBqX +XiJ,Bfx —fiAfx, 
Aly = -XBly -^I^B*y "MS?/- 

Applying any of the two parts of Theorem [3] yields an associated GEP corre- 
sponding to a matrix pencil method. As an additional permutation, the order 
of the two equations in a two-parameter eigenvalue problem does not influence 
the problem and can be interchanged to yield yet other variants. Hence, in 
total for the neutral single delay DDE we find 4 ■ 2 ■ 2 = 16 matrix pencil 
variants. (For DDEs with multiple delays there is one associated polynomial 
eigenvalue problem resulting in 8 variants.) 

The methods known in the literature correspond to the following: 

• and \S\ : fl26p and the A-part of Theorem [31 

• [T9j : (|T9l) and the //-part of Theorem [3l 

• [H, [I3j, and ^: (dnD and the A-part of Theorem O 

Finally, we stress that the above list, which contains many new variants, is 
more than just an theoretical encyclopedic description of all possible options. 
Depending on the given matrices, the structure and sparsity patterns of the 
Kronecker products may differ which may imply that for certain applications 
some methods may be more favorable than others. 
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6 Conclusions and outlook 

We have recognized new types of eigenvalue problems: quadratic and polyno- 
mial two-parameter eigenvalue problems. Using these problems as a unifying 
framework, we have derived associated (one-parameter) quadratic or polyno- 
mial eigenvalue problems that are at the heart of many matrix pencil methods 
that are used to analyze asymptotic stability of DDEs. This unifying way to 
derive the matrix pencils in the matrix pencil methods provides further under- 
standing of these methods and makes it easier to compare various approaches. 
Moreover, we have proposed several new variations on known matrix pencil 
methods. 

Furthermore, we expect that the recognized framework of quadratic and poly- 
nomial two-parameter eigenvalue problem may lead to a considerable amount 
of new research. First, we want to stress that it has been outside of the scope 
of this paper to study theoretical and practical properties of these new types 
of eigenvalue problems. There are many interesting aspects that need further 
investigation, such as how to carry over the concept of linearization (as is 
common practice for QEPs and PEPs) to these problems. 

Second, the matrix pencils constructed by Kronecker products that occur in 
the matrix pencil methods are of large dimension by nature, even for medium- 
sized problems, which may make efficient computation of eigenvalues and sta- 
bility of DDEs very challenging. We believe that the key to a successful com- 
putational approach lies in a direct attack of the polynomial two-parameter 
eigenvalue problem, instead of the corresponding matrix pencils, in the same 
spirit as, for instance, [12] for the linear two-parameter eigenproblem. We 
leave both of these topics for future work. Third, we note that an upcoming 
work by Muhic and Plestenjak will examine relations between quadratic two- 
parameter eigenvalue problems and singular linear two-parameter eigenvalue 
problems [23] . 

Acknowledgment: The authors would like to dedicate this work to Prof. Henk 
van der Vorst, who has always pointed us to new interesting and challenging 
problems and has provided inspirational ideas to deal with them. 
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